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Abstract 

Linear exact modeling is a problem coming from system identification: 
Given a set of observed trajectories, the goal is find a model (usually, a 
system of partial differential and/or difference equations) that explains 
the data as precisely as possible. The case of operators with constant 
coefficients is well studied and known in the systems theoretic literature, 
whereas the operators with varying coefficients were addressed only re- 
cently. This question can be tackled either using Grobner bases for mod- 
ules over Ore algebras or by following the ideas from differential algebra 
and computing in commutative rings. In this paper, we present algorith- 
mic methods to compute "most powerful unfalsified models" (MPUM) 
and their counterparts with variable coefficients (VMPUM) for polyno- 
mial and polynomial-exponential signals. We also study the structural 
properties of the resulting models, discuss computer algebraic techniques 
behind algorithms and provide several examples. 
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1 Introduction 



Linear exact modeling is a problem of system identification that leads to in- 
teresting algebraic questions. We start with some motivation from the systems 
theoretic point of view: The problem of linear exact modeling was formulated 
for one-dimensional behaviors in [T] , see also [9J [TU] . Starting with an observed 
set of polynomial-exponential signals, the aim is to find a linear differentiation- 
invariant model for these. Evidently, the whole signal set is a behavior that 
is not falsified by observation. But such a model has no significance. Making 
the behavior larger than necessary, the accuracy of the explanation decreases. 
So besides the condition that the desired model should be unfalsified, we are 
searching for the most powerful one. This means that the model does not admit 
more solutions than necessary. A model satisfying all conditions is abbreviatory 
called continuous MPUM (most powerful unfalsified model). 
In [TB] , the modeling was extended to multidimensional behaviors [3J [Tl] , and 
in [17 to the discrete framework, that is, instead of the requirement that the 
model should contain all derivatives of the signals, it is required that all shifts 
of the signals are contained. 

In other words, the problem is to find a homogeneous system of partial differ- 
ential equations with constant coefficients that is as restrictive as possible with 
the property of possessing the observed signals as solutions. 
In [T5] a different approach was introduced. There the goal is to find all partial 
differential equations with polynomial coefficients that are solved by the signals. 
Thus the new aspect of this approach is the choice of a different model class. 
Indeed the properties of the resulting model depend strongly on the model class. 
For instance, by the transition from the MPUM to the VMPUM, the time- 
invariance vanishes. In this paper, we continue this approach. But since the 
continuous case is not the only interesting one, we will consider a more general 
problem comprising both the continuous and the discrete situation. Later some 
special model classes will be discussed in more detail. 

Let us particularize our goal. Let K be a field and O be an operator algebra over 
K . Further let Ao be a function space over K possessing an O-module structure. 
A model or a so-called behavior B is the solution set of a homogeneous linear 
system, given by finitely many equations. These equations are defined in terms 
of the operator algebra O. Thus B is characterized by 

B = So\(0 lxr R) = {uj£ A% \R»u = 0}, where R G O rxm 

and • denotes the natural extension of the module action o • u) of o € O on 
oj G Ao to the matrix R G O rxra and the vector oj G A l . In most cases of 
interest, we have K C O and ok — ko for all k G K, o 6 O. Then B is a 
if-vector space, and thus the introduced model class is linear. Within such a 
model class we want to perform modeling now. Suppose we observe a set of 
signals f2 C Aq. The aim is to find a model Bn in the model class such that 

1. Bq is unfalsified by tt, i.e. fl C Bn- 

2. £>q is most powerful, i.e. for every behavior B with fl C £>, it follows that 
B n C B. 

If Bn is invariant under the action of O, that is, if we have for all o G O 

uj G B => omoj G B, 
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it is called most powerful unfalsified model, short MPUM of f2. Else, if Bq 
varies under it is called variant most powerful unfalsified model, short 
VMPUM of 0. We denote the VMPUM of fi by B%. 

The following example shows how the choice of the model class affects the model. 
Example 1.1 Consider the signal set consisting of a single signal 
fi = {u}, where uj(t) =t for all t€R. 

1. Let O = C[d] and Ao = C°°(JR, C), where d • / := f t . Using the commu- 
tative structure of the operator ring, the underlying system is invariant 
under differentiation: 

R • w = R(d • to) = (Rd) »w= (dR) •w = d(R •w)=0. 

Since we are searching for a differentiation-invariant model, we obtain 
that besides u, also its derivative, the constant function 1, belongs to £>n- 
Using that the model is C-linear, we get that 

B n = {to | 3a, b G C : Vt G R : w(t) = at + b}. 

An element w G C°°(K, C) is contained in Bn if and only if 

d 2 • to = 0, 

i.e. the MPUM is specified by a single ordinary differential equation with 
constant coefficients. 

2. Now let O — C[t]{d), where d • / := ^ and Ao is defined as above. 
We want to describe w as a solution of homogeneous ordinary differential 
equations with polynomial coefficients. The equations 

d 2 • to = and td • to — w = 

are satisfied by lj. We will see later that these two generate a kernel 
representation of the VMPUM of f2. The corresponding solution space 
equals 

B^ = {w | 3a £ C : Vt £ K : tw(t) = at}. 

Notice that this example demonstrates the variance under d, since we 
have d»u> £ Bq. Another property that should be pointed out is that the 
VMPUM yields a more precise description of fi than MPUM. 

2 Ore algebras 

The example above deals with continuous signals. But in applications, there 
are also discrete phenomena or combinations of discrete and continuous signals 
that are of great interest too. Many of the relevant operator algebras have the 
structure of an Ore algebra, as studied e.g. in [3J [21 H]. We give a definition 
that is motivated by g]. Moreover, this simplifies more general setup of [5]. 
Hence first consider skew polynomial rings, a generalization of polynomial rings 
to the noncommutative framework. 
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Definition/Remark 2.1 [T3] 

(1) Let A be a ring and a : A — > A be a ring endomorphism. 

(a) The map 5 : A — >• A is called a er-derivation if it is /^-linear and 
satisfies the skew Leibniz rule 

6(ab) = a(a)6(b) + S(a)b for all a, b e A. (I ) 

(b) For a u-derivation 5, the ring A [d; a, S] which consists of all poly- 
nomials in d with coefficients in A with the usual addition and a 
product defined by the commutation rule 

da = a(a)d + 5(a) for all a G A, 

is called a skew polynomial ring or an Ore extension of A with 
a and 6. 

If A is a domain and a is injective, the skew polynomial ring A [d; er, S] is 
a domain by degree arguments. Then the definition can be iterated to the 
so-called Ore algebras [I]. 

(2) Let A = K [ti, . . . , t n ]. An iterated skew polynomial ring 

O = K[ti, . . .,t n ][di;ai,6i] ■ ■ ■ [d s ;a s ,5 s ] 

is called a (polynomial) Ore algebra if the oVs and 6j's commute for 
1 < z , j < s , the cVs commute with dj's and further for all 1 < i < s the 
map er, : O — > O is an injective if-algebra endomorphism and Si : O — >• O 
is a (Ti-derivation satisfying (Ti(dj) — dj and Si(dj) = 0. 

Using multi-index notation, every element of an Ore algebra can be expressed 
into the normal form 

J2 Pad a = Yl P*di ai ■■■■■ 'V where Pa e A. 

For our issues the most interesting examples of Ore algebras are the following 
ones. 

Example 2.2 Let n = 1, thus A = K[t}. The algebras can be iterated to 
n e N. 

1. The first Weyl algebra is defined by Wi := A[d; idwi i ~§i] w hh the com- 
mutation rule dt = td + 1. 

2. The first difference algebra is defined by Si := A [A; cr, 6], where (ap)(t) = 
p(t + I) and 5(p) = cr(p) — p for all p e 5i. The commutation rule is 
At = tA + A + 1. 

3. The following Ore algebra is a combination of the first and second one. 
Define SWi := A[A;ai,Si] [d; (72,62], where 02 '•= idswn 62 '■= §1 and 
{(7ip)(t) = p(t + I), Si(p) = o\ 0) -P for allp G Then 9t = td + 1, 
At = tA + A + I and 9A = Ad. 



4 



4. Suppose q to be a parameter. The first continuous (^-difference algebra 

is denned by Q := A[d; a, 5], where a(p) = p(qt) and 8{p) = p(qt) — p{t). 
We obtain the commutation rule dt = qtd + (q — l)t. 



Lemma 2.3 Let A be a ring, and A[d;a, 5} be an Ore extension of A. For any 
a G A there exists an Ore extension A[A a ; a, 5'] with 5' (a) — a(a)a — aa + S(a), 
such that A[d; a, 5] = A[A a ;a, 5'] as rings. 

Proof For all a E A, the equality da — a(a)d + 5(a) holds. For a £ A define 
A Q := d — a. Then it obeys the relation A a a = a(a)A a + <j(a)a — aa + 5(a) = 
a(a)A a + 5'(a). The map 5' is linear and it is a cr-derivation since 

5'(ab) = a(a)a(b)a — a(a)ab + a(a)5(b) — a(a)ab — aab + 5(a)b = 

a(a)a(b)a — aab + a (a) 5(b) — 5(a)b = a(ab)a — aab + 5(ab). 

Define the ring homomorphism tp a : A[d;a, 5] —5- A[A a ; a, 5'], tp a (a) = a for all 
a G A, ip a (d) = A a = d — a. Then ip a is an isomorphism. □ 

Let O := A[di;ai,5\] ■ ■ ■ [d m ; a m , 5 m ] be an Ore algebra. With the action 

di*p := Si (p) and amp := a- p for all p G A and a e A 

the if-algebra A becomes an O-module. For this, we have to show that 

1. (oi • 02) mp = o\ m (02 m p) for all o\ , o 2 G O and p E A 

2. (01 + 02) mp = o\ m p + o 2 mp for all 0\ , o 2 G O and p G A 

3. o • (p + q) = o m p + o m q for all o G O and p, q G A. 

To show 1. it suffices to consider 01 = adi and 02 = bdj with a, b G A. Then 

(o 1 -o 2 )mp = (a(ui(b)di + 5i(b))dj)mp=(aui(b)didj+a5i(b)dj)mp 

= aui (b)5i (p)Sj (p) + a5 t (b)5j (p) — a5 t (bSj (p)) = adi • ( bdj • p ) 
= o 1 m(o 2 mp). 

The equality in 2. and 3. holds by similar arguments. 

Using this action, we can define the kernel of a linear operator / from the Ore 
algebra O over a ring A to be ker^ / := {a G A \ f ma = 0}, which is a K- vector 
space. 

Lemma 2.4 Let K be a field, A be a K -algebra, d be a K-linear operator, 
acting on A and B = A[d;a, 5] be the corresponding operator algebra (that is, 
for all a G A we have da = cr(a)d + 5(a)). Then the following holds: 
(1) \cy a d = A^-5 = {)^B = A[d; a, 0]. 

(ii) If kcvA d — A, then we have for A := d — 1: A[d; a, 0] is isomorphic as K- 
algebra to operator algebra A[A; a, 5'] with 5' := a — 1. Moreover, ker^ A={c£ 
A I a (a) = a} = const^ A Q A with the equality if and only if A is invariant 
under a, what is the case if a = 1a- 
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Remark 2.5 Using Lemmata 12.31 and 12.41 we pass to the new setting of opera- 
tors, which action • is nontrivial on A. We call such an operator nontrivial and 
from now on, we work with such operators only. 

Example 2.6 Consider the two most important operator algebras, built from 
operators having zero kernels. The first forward shift algebra is defined by 
K[t] [s;a,0] with (af)(t) = f(t + 1) for all / £ K[t}. The commutation rule is 
st = ts + s. There is a natural operator associated to s, namely the difference 
operator A = s — 1, already defined in !2.2[ obeying the relation At = tA+A+1. 
Applying Lemma 12.41 we see by degree argument, that ker A = K and the two 
algebras are isomorphic both as Ore extensions and iC-algebras. 
Let q be transcendental over K. Then the first ^-commutative algebra (or 
Manin's quantum plane) is defined as K q [x, y] := K(q)[x] [d; a, 0] with (af)(x) = 
f(qx) for / £ K[t]. Again, there is a natural q-diffcrcncc operator A q := d — 1 
and the corresponding operator algebra has been already described in 12.21 as 
the first continuous g-difference algebra. Its commutation rule reads as dt = 
qtd+(q-l)t. 

For oi, ... ,Ok £ O n , we denote by o(°i, ■ ■ ■ °k) the left submodule of O™, gen- 
erated by o±, ... ,Ok- 

Theorem 2.7 Let O be an Ore A-algebra, built from operators d\, . . . ,d s which 
have non-zero kernels. Then there is an isomorphism of left O-modules 

0/o(9i, ■ ■ ■ ,d s ) = A. 

Proof There is a left O-module homomorphism 

<p : O A, a= a a d a i-f a • 1 

since ip(b - a) — (b ■ a) • 1 = b • ip(a). Due to Def. 12.11 ([T]) we have 5(1) = 
and thus a • 1 = ag. The kernel of if is given by the left ideal o{di, . . . ,d s ). 
Further, ip is clearly surjectivc. So the claim follows from the homomorphism 
theorem. □ 

Following Theorem 12. 71 every polynomial p £ A can be viewed as an element of 
the left O-module Ofo(di, . ■ ■ , d s ) by identifying p with p + o (di , . . . , d s ) =: [p]. 
Then the action of di is exactly the a^-derivation 5i, since 

di\p] - [d iP ] - [o-Mdi + 5i(p)] = Mp)} = [di.p]. 



Remark 2.8 Let p £ A and o £ O. Then there is the following equivalence 

o • p = if and only if o ■ p £ o(d\, . . . , d s ). 

Proof By Theorem l2.71 we have an O-module isomorphism A = Ofo(di, ... , d s ) 
given by 

A^0/ o (di,...,d s ), p^[p\. 

Since the O-module structure is respected, o • p maps to [o ■ p] and hence the 
claim follows. □ 
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Remark 12.81 gives the possibility to describe and to compute the annihilator of 
an element p G A. Consider the map 

K p : O -> 0/ o (di,...,d a ), o4o-[p], (2) 

which is clearly a left O-module homomorphism with the kernel 

ker(Kp) = Anno(p) := {o G O \ o*p = 0}, 

which is a left ideal in O. See Corollary 13.21 for its algorithmic computation. 
This construction lifts to the case of vectors. Suppose p — [pi, . . . ,Pm] T £ A m . 
An element of o 6 lxm naturally acts on p by 

I n 
i=l 

A subset B C A m is called invariant under G C lxm if and only if o »p = 
for all o 6 G and p € B. The set of elements under which p is invariant has an 
O-module structure and equals to the kernel of 

m 

k p : lxm -> 0/ o (di, . . .,9 S ), o = [oi, . . . ,o m ] ^ 0i • [p l ]. 

Moreover the following isomorphism holds 

lxm /ker(K p ) = o(pi, ■ ■ ■ ,Pm)/o{Pi, ■ ■ ■ ,Pm) ri (di,...,d s ). 

The image of k p equals (o(pi, ■ ■ • ,Pm) + o(^i, ■ ■ -,d s ))/o{di, ■ ■ -,d s ). This is 
isomorphic to o(pi> • • • ,Pm) / o(pi, ■ ■ ■ iPm) ^o(9i, ■ ■ ■ , 9 S )- So the claim follows, 
since k p is a homomorphism. 

Remark 2.9 If O is Noethcrian (see [3]), then the left submodule ker(ft p ) C 
lxm is finitely generated. 

For a polynomial m-tuple p G A m , we consider 

Ann (p) = {oeO|o.p = 0} = {oeO|otp 1 = 0Vj} = f]Ann (p t ), 

which is a left ideal in O. As we see immediately, Anno(p) lxm is a (usually 
strict) submodule of ker(K p ) and hence, the latter typically has more interest- 
ing structure, see Example 16.31 It is always possible to recover Anno(p) from 
ker(Kp). In our opinion, using ker(/-c p ) is more natural in the context of vectors 
of signals. 

3 Algorithmic computations 

For the concrete calculations used in this article, we need algorithms for the 
following computational tasks over (polynomial) Ore algebras: 

1. syzygy module of a tuple of vectors 
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2. elimination of module components from a submodule of a free module 

3. annihilator ideal of an element in a finitely presented module 

4. kernel of a homomorpliism of modules 

5. intersection of a finite number of submodules of a free module. 

Let O be a Noetherian Ore algebra. Moreover, let M be a finitely presented 
left O-module, that is, there exists a matrix P £ £) mxn such that there is the 
following exact sequence of left O-modulcs: 

Qlxm ^ Qlxn ^ M ^ o 

Recall that for a tuple F = (fi, . ■ ■ , f s ), fi C lxn , the set LeftSyz(F) := 
{[at, . . . , a s ] £ lxs | J2i a ifi — 0} carries the structure of a left O-module and 
is called the left syzygy module of F. Since O is Noetherian, LeftSyz(F) is 
finitely generated. Computation of syzygies over Noetherian Ore algebras can 
be accomplished with several algorithms and requires Grobner basis techniques; 
see [5] for Ore algebras and [7] for the commutative case. 

£ 

Let {a} be the canonical basis of the free module lxl = @Oej. 



Proposition 3.1 1. "Elimination of module components" . 

Let S C lxi be a submodule. Moreover, let <o be a monomial ordering 
on O and < m = (c, <o) be a position over term monomial module ordering 
on the free module lxe , defined as follows. The components are ordered 
in a descending way e\ > ■ ■ ■ > eg and for any monomials 0\, 02 £ O 

0\&i < m Oitj 44> j < i or (j — i and 0\ <o 02)- 

Let G be a Grobner basis of S with respect to < m . Then V 1 < k < I 

I l 
Gn ffi Oei is a Grobner basis of S (~l © Oei . 

i—k i—k 

2. "Kernel of a homomorphism of modules". 

Consider an O-module homomorphism lxs % lxn /O lxm P, &i i_ » [^], 
where £ lxn . Let Pi be the i-th row of the matrix P. Then 

s 

ker^ = LeftSyz( P X) P m ) )n0Oe,. 

i=l 

e 

Proof 1. Define W = ^Oe^. Since G is a Grobner basis of S, for any 

i—k 

s £ S there exists g £ G such that lm(g) divides lm(s). If s £ SPiW, then 
lm(<?) £ W and hence, by definition of < m , we have g £ W and g £ GtlW. 
So, G n W is a Grobner basis of S D W. 
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2. We have 

s m 

[h,...,b s ] ekerip 3a k e O : ^ b^i + a k Pk = 

i=l fe=l 

s 

[6i,...,6j GLeftSyz( . . . , P x , . . . , P m ) )n0Oe,. 

»=i □ 



Corollary 3.2 1. "Annihilator of an element in a module" . 

Let M = lxn /0 lxm P and letP u ...,P m be the rows of P. Moreover, let 
v € lxn . Then the left ideal Ann£ r (v) := {a e O \ a[v] = e M } C O 
can be computed as 

Ann^(u) = ker(0 '4 1 M) = LeftSyz( («, Pi, ... , P m ) ) n Oei. 

"Intersection of finitely many submodules" . 

Let N±, . . . , iV m C lxr be submodules. Then 

m 

f]Ni = ker(0 lxr (0 lxr /JVi)©. ■ •©(O 1 ^™), e* ^ (N, . . . , [ej)). 



Remark 3.3 For an O-module homomorphism lx 70 lxr Q^0 lx, 70 1><m P, 
its kernel is the image of ker-0 (as in Theorem 13. ip under the natural projec- 
tion lxs — > lxs /0 lxr Q. A left Grobner basis can be obtained by reducing a 
left Grobner basis of kert/j + lxr Q with a left Grobner basis of lxr Q, see 

Note that in practical computations, elimination of module components is usu- 
ally not complicated. This stands in distinct contrast with the elimination of 
algebra variables, which is often very hard to achieve. The algorithms used 
in this article involve only the elimination of module components and thus are 
feasible in practice. 

The algorithms we have discussed are implemented in computer algebra systems 
like e.g. Singular:: Plural [5] or Maple [510] with the package OreMod- 
ULES. More background on these algorithms can be found in e.g. [8], [11] . 

In particular, a set of generators of ker(/c p ) from the previous section can be 
calculated explicitly. 

By Proposition I3.1f 2) ker(« p ) is obtained via the kernel of a module homo- 
morphism, that is, by one Grobner basis computation with respect to module 
monomial ordering eliminating components. The monomial part of this ordering 
can be chosen arbitrarily to be e. g. a fast one. 
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4 Application to linear exact modeling 



We will now use the results from above to define an unfalsified and most powerful 
model over an Ore algebra. 

Assumptions and notations: Suppose O to be a Noetherian Ore algebra 
with the additional property that di acts nontrivially on A for all 1 < i < s. 
Recall that Ao denotes a function space over K possessing an O-module struc- 
ture. Suppose further that A C Ao- 



Remark 4.1 [ [T3], Theorem 1.2.9.] Since A is Noetherian, O is Noetherian if 
tji is an automorphism for all 1 < i < s on A. 

Thus all Ore algebras considered in Example 12.21 are Noetherian. 

Suppose that the cr-derivation acts nontrivially on A. Then the corresponding 

most powerful unfalsified model varies under O. 

Starting with a single signal p 6 A m , we want to find the VMPUM of p, that is 
a behavior, invariant under some finitely generated submodule of lxm . 



Theorem 4.2 Let p 6 A m be given 
ker(Kp) = o(ki, ■ ■ ■ , K) and let R £ O rxm 
Then the VMPUM of {p} is given by 



Consider the map n p from 0). Let 
be a matrix whose i-th row equals ki. 



B\ p} = {g e A"S 



R»g = 0}. 



Proof By the definition of R and Remark 12.81 it is clear that {p} C BY p -,. 
It remains to show that ${ p } is most powerful. Suppose there exists another 
behavior £>' unfalsified by p. The behavior £>' possesses a kernel representation 
R' e O r ' xm . By the definition of i?, there exists a matrix X G O r xr such that 
R' = XR. But since (X ■ R) • p = X • (R • p), it follows that B\ p} Cg'. n 

Example 4.3 Let us consider a more interesting example than Example 11.11 
with respect to our favorite algebras from Example 12.21 Let ft — {u>} consists 
of the cuspidal cubic 

w(*i,*2) = 4-4- 

Let us denote by Ao — C[[fi,t2]] the ring of formal power series and consider 
the VMPUM BY, = {/ € Ao \ #vmpum • / = 0} of O with respect to several 
operator algebras O. 



1. Suppose O to be the second Weyl algebra (see Example 
using Singular we obtain: 



di 
did 2 



Then by 



R 



VMPUM 



df + 3<9f 
tidl - d 2 

t 2 dl + Mid 2 

2i 1( 9i + 3t 2 5 2 - 



Now let us determine BY u x to see how precise the description given by the 
VMPUM is. Let / e Ao- 
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(a) <9f • / = => f = Co + ci< 2 + c 2 t\, where a G C[[ti}]. 

(b) dxd 2 »f = dx»c 1 +2t 2 di»c 2 = => 9i »ci = A d x »c 2 = 

ci,c 2 e C. 

#i • c + 6c 2 = 



(c) (a 3 



1 - 



/ = o 



C() 



-c 2 if + rf 2 *i 



do, where di € 



(d) (t 2 a 2 2 - d 2 ) . / = o 



ci =0. 
=► d 2 = 0. 
= =>■ -4di£i 



(e) (3t 1 9 2 +t 2 9 1 2 )./ = 

(f) (2ti9i + Zt 2 d 2 - 6) • / = =>• -4diti-6d =0 d a = = do- 
Hence, we obtain that / = c(t\ — t^), thus 

#M - W*? - $ I c 6 C}. 

With respect to the requirement of being most powerful and linear, the 
VMPUM is as significant as possible. We observe that the VMPUM of a 
single non-zero signal has C-dimension one. Actually, this holds in general, 
as will be shown in Theorem 15.21 



2. Suppose O to be the second difference algebra see Example 
by using Singular we obtain: 



Then 



R 



VMPUM 



Ai 



AiA 2 



A? 



3A 2 



2i 2 A 2 -f aa 2 



8Af 



21A| 



A 2 - 2A 2 
2i 2 A2 + Af + 6iiA 2 + 6A 2 
+ 24t 1 A 1 + 36i 2 A 2 - 24Ai - 



18A 2 - 72 



Similar arguments as above lead us to 

= W*i - *a) I c G C}. 

3. Suppose O to be the second SW algebra see Example 1 2. 2 1 Then by using 
Singular we obtain: 

Ai 
AiA 2 
A 3 + 3A 2 



R 



VMPUM 



2<9 2 + A 2 2 - 2A 2 



2<9i + A( - 2Ai 



2A 2 



8Af 



2t 2 A\ + A| - 2A 2 
2t 2 Af + A2 + 6tiA 2 +6A 2 
21 A?, + 24*i A x + 36i 2 A 2 - 24Ai - 18A 2 - 72 



Note that generators in the output depend on the monomial ordering of 
the operators. In this example Ai j2 were chosen to be greater that di, 2 . 
Taking a reverse ordering produces different (but equivalent) answer. 

Comparing this matrix with the matrix above, we see that the rows of 
the matrix belonging to the difference case appear also here. We conclude 
that 



{c{t\ -4)\ce C}. 



n 



Thus, taking SW as operator algebra, we have got more equations than 
with the difference algebra. However, we have obtained very interesting 
mixed differential-difference equations, which show the interplay of two 
different operator settings. 



4. The second y-difference algebra see Example 12.21 

-RvMPUM = 



dl + (-q 2 + 1)52 
(-q - 1)9! + (-q 2 - q - l)8 2 + (g 4 + q 3 - q - 1) 

i{d 2 - tld 2 + (g 2 - \)t\ 



(a) The first equation yields: 

E - i)H\4 + (-q 2 + 1)^2 mi = o 

(q j - I) 2 + {-q 2 + -1) = 
& j = 0Vj = 2. 

(b) Now consider the second equation. 

i. Suppose j ' = 2, then 

(-9 - 1)5]^^ - 1)4*2 + (-1 2 - 1 - l)5>i(<? 2 - 



(g 4 + 5 3 -9-1)^^-4*1=0 



& (-q - l)(q l - 1) + (-q 2 - g - l)(q 2 - 1) + {q i + q 3 - q - 1) = 

i = 0. 

ii. Suppose j = 0, then 

(- q ^o(g % - 1)4 + (? 4 + q 3 - q - 1) E ^4 = o 

i = 3. 

Thus / = e 30 i? + 002*2- 
(c) Applying the last equation, we get 

t\c 2{q 2 - - <2C0 2 (g 1 - 1)*2 + - 1)*2(C304 + C 02 ^) = 

tl4(q 2 - l)(c 30 + c 02 ) =0 ^ c 30 = -c 02 . 
Thus we obtain once more = {c(tf — t 2 ) \ c E C}. 



Remark 4.4 As we have seen in the previous example, the number of equations 
giving the VMPUM depends strongly on the underlying Ore algebra. In all 
cases, with Grobner bases we get more equations than it might be actually 
necessary. However, it is possible to compute a smaller generating set, which is 
usually not a Grobner basis. Namely, one computes a left syzygy module of a 
given system and almost directly deduces a smaller generating set from it. As 
an example, we show that only 3 of 6 equations from the first example of 14.31 
generate the whole ideal, namely 8182, df + 3d 2 , 2t\d\ + 3t 2 d2 — 6. Analogous 
smaller generating sets can be obtained for other examples. 
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Theorem 14.21 can be generalized to a set of several signals directly. A kernel 
representation of the VMPUM of f2 = {oj\, . . . ,ujn} is determined by stacking 
a set of generators of 

N 

P| ker(K Wi ) 

i=l 

row-wise into a matrix R. 

Theorem 4.5 Using the notation from above, the VMPUM of fl equals 

BY 2 ={geA"S\R.g = 0}. 

Proof By the definition of R, it is clear that fl C B^. Also the property of 
being most powerful follows by the same arguments as used in the proof of 
Theorem g^] □ 

Example 4.6 Suppose O to be the first Weyl algebra and Ao = C°°(K,C). 
Consider the signal set f2 = {t, v$t — v\t 2 }, where vq, Vi G C \ {0}. The second 
trajectory will appear in Example 15.41 again. Since 

ker(« t ) n keT(n Vot _ Vlt 2) 

= Wl (td - 1, d 2 ) n Wl (~v 2 d 2 + {Av\t - 2v v 1 )d - 8v 2 l7 d 3 ) 
= wAt 2 d 2 ~2td + 2,d 3 ), 

the VMPUM of fi is given by = {c x t + c 2 t 2 \ ci,c 2 G C}. The intersection 
of submodules of a free module over a Noetherian Ore algebra can be computed 
as in Corollary 13.21 for instance with the system SINGULAR: : Plural [B]. 

5 VMPUM by using the polynomial Weyl alge- 
bra 

In this section, we suppose O to be the n-th Weyl algebra 

O = W„ := C[ti, . . . ,i„][<9i;id w „, ^-] • • • [d n ; idw„, 

Thus for p G C[ti, . . . ,t n ], we obtain di • p := J^. Further suppose Ao to be 
C°°(R n ,C), the space of smooth functions. Identifying a polynomial with the 
corresponding polynomial function, we obtain A C Ao- 

In this context, the VMPUM was already introduced in [TS]. Here, we will recall 
some results and additionally point out a new interesting property. 

5.1 C-dimension 

A known result is that the VMPUM is a finite-dimensional vector space over C, 
since it is contained in the corresponding MPUM [15] . In some cases, we can 
determine the dimension more precisely. We claim that the VMPUM of a single 
non-zero signal has C-dimension one. 



13 



Suppose p G A m . Every polynomial pi can be written as J2 k =i c ikt l3ik , where 
Cjfc e C for all i, k. Let £i := {/3a, • • • , Pihi} C Nq denote the set of all exponent 
multi-indices occurring in pi and let 

dij := max {(/3 ik )j \ /3 lk £ £*} (3) 
i<fc</ii 

be the highest degree in t,- of pi . Recall that by e; we denote the t-th canonical 
generator of the free module A m . The set 

£ Pz = { a E Mg- 1 ctj < dij + 1 for 1 < j < n} 

is finite, that is, £ Pi — {an, . . . , a^^ ^ }. Define for p 6 A m 

uer p - ^ qii ,...,p m , gamlm }■ 

Let Syz(Der p ) denote the module of polynomial syzygies. Define for the matrix 
M = [lei, d Qll ei, d a »iei, le m , d a ^e m , d amlm e m \^ the A- 

module homomorphism 

$ p : Syz(Der p ) -> ker(K p ), (gi, . . . , ^ (<?i, • ■ • , <7i) • M, 
which is clearly injective. 

Lemma 5.1 w n (Im($ p )) = ker(K p ). 

Proof Evidently w„(I m ($ p )) C kcr(K p ). Now suppose that a S ker(K p ). Since 
every element in W„ can be written in normal form, we obtain 

a • p = ^2 ak • "Pk = ^2 £ Ck i tPk] d1 " 3 ) ' Pk 

k k j 

k j 

Let us split the element a in a z and a nz such that a = a z + a nz and (a z )k 
consists of the parts of a k where d lkj • pk is zero. 

By the choice of dij, the set {dj** 3 * \ l < j < n, 1 < i < m} generates the 
set of d 1 with the property that there exists 1 < i < m such that d 1 • pi — 0. 
Then a z is contained w n {® * 1 1 J < n > 1 < * < rn)- But by the choice 

of Der p , the element a z is in the image of $ p . Suppose d lki »pk ^ 0, then 7^ is 
equal or smaller than (dki, ■ ■ ■ , dkn) in each component and again by the choice 
of Der p , the element a nz is contained in the image of $ p . Thus it follows that 
a e Im($ p ). □ 



Theorem 5.2 The VMPUM ofp 7^ is a one- dimensional vector space over C. 

Proof We use the notation of the Lemma 15.11 which reduces to commutative 
calculations. It is easy to see that the equivalence 

S '(fi,d a ".f 1 ,...,f m ,...,d a ^.f m f = ^ $ p ( s )./ = (4) 
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holds for every s G Syz(Der p ). Now let us discuss the left hand side. More 
precisely, let us consider the solution space Sol(Syz(Der p )) in Ao belonging to 
Syz(Der p ). Since Der p contains all non-zero derivatives of pi for all i, there 
exists a non-zero constant C 3 k 6 Der p . We can suppose k — 1 and without 
loss of generality let d amhm • p m — 1. Then 

(-l,0,...,0,Pi), (0, ■ ■ ■ , 0, -1, ) G Syz(Der p ) 

and thus 

Sol(Syz(Der p )) = {c • (p x , . . . , ^J" ) | c G C}. (5) 

Now suppose / G VMPUMp. From Lemma [5.11 together with ((4]) and ([5]) we 
deduce the claim. □ 



Remark 5.3 In the case of a single non-zero signal, the VMPUM gives the 
most precise description one can get for a linear system. 

Example 5.4 Consider the trajectory uj(t) = vot — vit 2 , where Vq, v\ G C\{0}. 
Then the MPUM of lj is given by 

B {u} = {a(v Q t - Vl t 2 ) + p(v - 2 Vl t) + 7(-2wi) | a^jel} 
= {at 2 + bt + c | a, b, c G R} 
= {w G ylo I <9 3 •«) = 0}. 

Thus there are three free parameters to choose. The VMPUM of cj is given by 



vld 2 + {Avft - 2v v 1 )d ■ 



3 



= 0} 



= {w G Ao 
= {c (v t - v x t z ) | eel}, 
that is, two degrees of freedom vanish when we consider the time- variant model. 



5.2 Structural properties 

Let us discuss some structural properties of the W„-module ker(K p ). Since every 
element of W„ can be transformed into normal form, the degree of an element 
a G W ra can be introduced as 

n 

deg(a) := max{ an + pi \ a = a a ^t a d^ , a a ^ G C}. 

i=l a,/3£N™ 

Then J^(W n ) := {a G W„ | deg(a) < i} induces a filtration on W n . The corre- 
sponding associated graded ring Gr(W„) is isomorphic to C[ii, . . . , t n ,d\, . . . , d n ] 
as a graded C-algebra. For every finitely generated W„-module M, we define 
the Hilbert polynomial HP]JJ" := HP^.^^. The dimension of M is defined 
as dim\v„(A/) := deg(HPM) + 1. Furthermore M is called holonomic if it has 
dimension n. A holonomic module is of minimal dimension, since the dimension 
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of W„-modules is bounded below by n and bounded above by 2n. Holonomic 
W„-modules are additionally cyclic and torsion modules. For details see [5]. 
As usual, we write A = C[ti, . . . , t n ] c W„. 



Theorem 5.5 There is an isomorphism ofW n -modules W n xm / ker(n p ) = A 
In particular, W n xm / ker(/-c p ) is simple holonomic W n -module. 

Proof Since n p is a homomorphism of W„-modules, we get 

W n xm / ker( Kp ) = lm( Kp ) C W„ / w „ (ft, ... , 9 n ) S A. 

Thus W^ xm / ker(Kp) is isomorphic to a submodule of A Due to the fact that A 
is a simple holonomic W„-module and W n xm /ker(K p ) ^ the claim follows. □ 



Corollary 5.6 Since W n xm /ker(K p ) is holonomic, there exists a left ideal L p , 
depending on p, such that W n xm /ker(K p ) is isomorphic to the cyclic left Wn- 
module W n / L p . 

An algorithm, using Grobner bases, to compute a generator of a holonomic 
module is given in [12] . On the other hand, since W* xm /ker(K p ) is simple 
holonomic module, any non-zero element can be taken as a generator for a 
cyclic presentation. 

Example 5.7 Suppose uj — [ci, C2, c 3 ] T for ci,C2,c 3 6M \ {0}. Then 
ker(K w ) = Wl ([0,c 3 ,-c 2 ], [c 3 ,0,-ci], [0,0,5]). 



Since 



C 3 — C2 

c 3 -ci 
5 



l/c 3 ci/c 3 
l/c 3 c 2 /c 3 
1 





" 1 





" 







1 













d 



we obtain Wf /kev( Ku ) = W? /ker(/c u )C ^ Wi / Wl (9) = C[i]. 



5.3 VMPUM of polynomial-exponential signals 

In this section, we extend the signal space that should be modeled. The goal is 
to compute the VMPUM of 

p = [pi exp A i , p m exp Am ] T , (6) 

where for all 1 < i < to, we have pi G A, A 1 € C™ and 

exp A := exp(Ai<i + • • • + X n t n ) for A e C™. 

By the action dj • exp A = Aj exp A for all 1 < j < n the space of polynomial- 
exponential functions becomes a W„-module. 
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Consider the scalar setting first, that is, m = 1. Define for A G C" the W ra - 
homomorphism 

o-a:W„->W„, di\^{di-\i), ti- 
lt is easy to see that er\ is a W n -automorphism. We claim that for a G W„ and 

feA ' 

a • p = if and only if a\(a) • (pexp A ) = 0. (7) 
For the proof suppose a = at c ' i d^ i . 

Using the identity {di — A») • (pexp A ) = {di • p) exp A , the claim follows by 
a A (a).(pexp A ) = ^ Ci ^((ai-Ai) Al ---(9 n -A„)^»).(pe XPA ) 

= ^at^ddf' 1 • • • <9*" ) • p) exp A = (a»p)exp A . 

Extending the dimension, there are two special cases requiring attention. First 
suppose A 1 , . . . , A m to be equal, that is, p = \p\, . . . ,p m ] T exp A , where A := A 1 . 
Then claim ([7]) can be generalized directly and it follows that 

m 

2J cti • (pi exp A ) = if and only if [a x , . . . , a m ] G cr A (ker(K p )). (8) 

i=l 

Assume now that A , . . . , A m are pairwise different. Then 

m m 

^ Qj • fa- exp AJ ) = if and only if [ai, . . . , a m ] G @ <t A j (ker(K Pj .) ). (9) 

Since exp A i , . . . , exp Am are algebraically independent over A, the claim follows 
from 

m 

E a i * ex P^ ) = 

3=1 

m I h ] \ 

E + Ai) (/3 - h ••• (d n + .pi exp AJ = 

3=1 \i=l J 

f7 Aj (dj) G ker(K Pj . ) for all 1 < j < m. 
Recapitulating we get: 

Theorem 5.8 Let f be of the form 0). Further let 

Ki := {j | X J ; = A*} ={k a ,...,k ih } 

and let I be chosen minimal such that we have a disjoint union KiO . . . UKi = 
{fen, . . . , k\h 1 , . . . , ka, . . . , hh t } — {!,..., to}. Further define the vector hi := 
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[fkn , ■ ■ ■ > fku .] T and Hi :— cr^i (ker(K/ li )). Let denote the kij-th canonical 
generator of W^ xm for 1 < i < I and 1 < I < hi. Defining for 1 < i < I, 

hi 

4>i : Hi -> W n , [a%, . . . , a hi ] h-> Oje fci3 , 
f/ie VMPUM o// is .gwen &y 0* =1 &(ffi). 

Proof After choosing a suitable projection, the claim follows by ([7]) and ([8"]).n 

6 VMPUM via the polynomial difference alge- 
bra 

Suppose that if = oo. Recall the definition of the n-th difference algebra: 
S n := K[ti, . . . ,t n }[Ai,ai,Si] ■ ■ ■ [A n ;a n ,S n ]. 

For p e K\t\, ... ,t n ], we have 

A, »p = 5i(p) = ai(p) ~p = p(t + a) -p(t). 

Further suppose that Ao = K N ° ■ Identifying a polynomial with the correspond- 
ing polynomial function, we obtain A C Ao- 

Similar to the continuous case, the kernel of be computed in a completely 

commutative framework. For this we choose a special representation of the 
polynomials that is adapted to the action of A, see [T7j. For t G Nq and 
v = (y\, . . . , v n ), we consider the binomial functions 

where (p) = 1 for all i. Then v\p„ = t\ ■ ■ ■ (ti — vi + 1) • • -t n ■ ■ ■ (t n — v n + 1) 
and moreover, each element p £ A m can be written as 

P= X! CuPv ( 10 ) 

for g G Nq, some suitable coefficient vectors c v G K rn and < cw denoting the 
component-wise order on Nq, that is, Vi < Qi for all 1 < i < n. Let us describe 
how to find this representation. We restrict to the scalar and one-dimensional 
case, where m — n — 1. The general case can be treated similarly. For p G A = 
K[t] we show how to find the introduced representation. Usually, a polynomial 
p is given in the form 

p(x) = d v t v + d„_ir _1 H h dit + d , where d. L G K. 

To write p in the form (|10[) , the occurring coefficients c v have to be determined. 
We will show how this can be done for a monomial d v t v . Since v\p v — t ■ (t — 
1) • • • (t — v + 1), we define 

9 (v) := t ■ (t - 1) • ■ • (t - v + 1) = t" + _1 + ' ' ' + 9[ U) t. 

First, the coefficients gi^ will be determined for 1 < v < v by using the fact 
that gM =g( u ~ v > -(t-v + 1). 



IS 



1. Determine g[^: 

The polynomial is a multiple of t. Recursively, one gets that 



1 for v = 1 

c-i)"- 1 rifc=i * fon/>i. 



2. Determine 



i/-2 



Using = • (ii/ + 1), we get 

= ffl"- 1) -("-i)-^- 1) = (-ir a n*-( i '- 1 )a 

*;=i 

Since = 1, we get a recursive formula. 

3. Determine for j < j/: 

A similar consideration as in the previous point yields 

ff] =s}-i - 1) -.9} • 

Finally, we observe 

<U B = d v (g{v)-g ( :l 1 -g{v-l)-{g ( :l 2 -g ( :l 1 .g ( :- 2 1) )g{v-2)--.-) 

I V-l \ / t)-l 

= d„ I + ^ fc„(i) - i) J = d„ I v\p v + ^ fc„(i) • (v - i)\ -p., 



where 

fc„(l) := -5^1, and k v (l) 



0, if I > v. 



Consider for example p(t) = t 3 + t 2 + 1. The bounding value g equals three, 
by using 



3 


pf 




Pf 




Mi) 


1 


1 








3 


1 


2 


-1 


1 





1 





3 


2 


-3 


1 









we finally get 

t 3 -- 

t 2 -- 

1 = 

=► p(t) = 



p(3) + fc 3 (l) • p(2) + fe 3 (2) • p(l) = 6 • p 3 + 3 • 2 ■ p 2 + 1 • Pl 

p(2) + fc 2 (l)-p(l) = 2-p 2 + l-p 1 

Po 

6 • p 3 + 8 • p 2 + 2 • pi + po 



In the following we show the advantage of this notation. Since 



{kp Vi ){ti) = 



ti + 1 

Vi 



! 2 — 

if Vi = 

U-i) if"i>l 
if = 0, 



19 



one gets, by using the fact that ^p v — S^p Ul ■ ■ ■ S^p Vn , the equality 







if P <c™ v 
otherwise. 



(11) 



Remark 6.1 Let p = (pi, . . . ,p m ) T 6 A m with Pi{x) — ad i it ,ldi + 
using multi-index notation. Dchnc 

Qi = max{(ui, ...,«„) e Ng = (/iki)j for 1 < fc < . 

Then the bounding multi-index g belonging to the binomial representation (|10[) 
is given by 

g = max{(vi, . ..,«„) u< = (^)i for 1 < j < m} . 

cw 

From now on suppose that p = X^eNj, v< cwB c vPv 

Remark 6.2 Connecting remark l6~Tl and (jlll) we get that S^p — for all \i 
with fii > Qi for at least one 1 < i < n. Now consider the finitely generated left 
A-module generated by 

Shiftp = a(S^P I p <cw £>)• 

The corresponding syzygy module Syz(Shift p ) is hnitely generated too, since 
A is a Noetherian ring. Analogously to the continuous case, we can give an 
A-module homomorphism from Syz(Shift p ) to ker(K p ), such that the image of 
Si, . . . ,Sd under this map generates ker(fc p ), that is, ker(«; p ) is finitely generated 
as an ^4-module. This implies that ker(/i p ) is hnitely generated as an iS„-module. 

Example 6.3 Let p = [t 3 ,t] T . Then the continuous VMPUM is the same as 
the discrete VMPUM, that is, equal to {c[i 3 ,i] T | c £ K}. Direct computation 
over Si yields 

ker Sl ( Kp )= 5l ([0,A 2 ], [0,*A-1], [1,-t 2 ]} 
and this means that 

" A 2 

tA-1 

1 -i 2 

is a kernel representation of the VMPUM of p. Note that over A\ , we have 

kev Al (K p ) = Al ([0,d% [0,td-l], [1,-t 2 ]}. 

Alternatively, we can compute ker^ (k p ) in the commutative framework, using 
the analogue of "difference algebra" approach. At first, we observe that 

Shift [t 3 )t p = K[t] (t 3 , 3t 2 + 3t + 1, 6t + 6, 6, t, 1) 



so 



- o - 




■ o - 


































1 









1 


1 







1 














. -f _ 




. -6 _ 




. -6t - 6 _ 





Syz(Shift [t3jt] r) = K [t]( 



Finally we get that ker (k p ) = Sl ([0, -tA+1], [A 3 ,-6A], [A 2 , (— 6i— 6)A], [A,(-3t 2 - 
3t-l)A], [1,-t 3 A], [A 4 ,0], [0,A 2 ])= 5l ([0,A 2 ], [0,tA-l], [1,-t 2 ]). 








" 1 " 


1 































- 3t - 1 _ 




-t 3 
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6.1 VMPUM of polynomial-exponential signals 

For A = (Ai, . . . , A n ) € K n the discrete exponential function is given by 

First suppose that m = 1, that is, we want to construct the VMPUM of a scalar 
polynomial exponential trajectory of the form pexp A , where p 6 A. Without 
loss of generality, we can assume Ai ^ for all 1 < i < n, since otherwise if 
Aj = 

pexp A (i) = j ° [J , where 

g . pjn-i t ^ (pexp A )(ti,...,ti_i,0,t i+ i,...,t n ). 

Consider the automorphism of S n 



ti I ^ tj 

A 4 ^ ^-(Ai-Ai + l). 



XA : $n ~> $n , 

Since the equality 

Xa(A;) • (pexp A ) = -^-(Aj - Ai + 1) • (pexp A ) 

= — (Ai exp A <Ti(p) - _pexp A -A t pexp A +pexp A ) 

= — (A i exp A (cr. i (p) -p)) = exp A A t »p 
M 

holds, we obtain the identity 

Xa(A?) • (pexp A ) = exp A Aj* • p 

that finally extends to 

Xa(A^) • Cpexp A ) = X ^(A) • (pexp A ) = exp A A^ .p. (12) 

Now using (fl~2j) we can deduce for a = J2i=i a ;A Qi G S n the equivalence 

a»p = & xx{a) • (exp A _p) = 0, (13) 

since 

h h 

X\(a) • (exp A p) = ^ajXA(A Ql ) • (exp A p) = ^a;(A Ql »p) exp A 

i=l i=l 
ft 

= cxp A 2J a« ( A Qi • p) = cxp A a • p. 
i=i 

Summarizing, we obtain 

Theorem 6.4 Lei i? € S l n xl be a kernel representation matrix of the VMPUM 
of p. Then the kernel representation matrix o/pexp A is given by (x\(Ri))i- 
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Now consider 



Pi exp A( i) 



p m exp A( „ 



where A (l) G K n \ {0} , and Pi e A (14) 



and suppose A (1 \ . . . , A (m) to be pairwise different and without loss of generality 
A^ l) ^ for all 1 < i < m and 1 < j < n. Then 

m m 

^ • (?>j exp A(j) ) = if and only if [ai, . . . , a m ] G (J) Xa«) ( ker(K Pj ) ), 

(15) 

which follows by 

rn ml hj 

a • (pj exp A(j ) ) = <S> £ £ Cjit^A^* • (p rf exp A0) ) J = Q 

m 

( X AU) ( a J ) ' ' ) ex PA<^> = Xxw ( a 3 ) e ker( Kpj ) for all 1 < j < 

i=i 

Choosing a suitable projection, we obtain by (fT5|) and (|15l) 



Theorem 6.5 Let p &e of the form (T^. Further let K l := {j | A J = A 1 } = 

{fcji, . . . , fey, } and / chosen minimal such that the disjoint union KiU . . . UKi = 
{fcn, . . . , kih±t • ■ • j ■ • ■ ? ^//iil — {1) • ■ • i m }- Further define the vector hi := 
[fk tl , ■ ■ ■ , fk a ] T an d Hi := x A (i> (ker(K^J). Let e^. denote the kij-th standard 
generator of S^*™ 1 for 1 < i < I and 1 < I < hi. Defining for 1 < i < I 

hi 

(pi : Hi ->• S n , [oi, . . . , a hi ] h-> ^ a j e k tj 

3 = 1 

the VMPUM ofp is given by 0' =1 

Conclusion 

Generalizing ideas from systems theory, we have defined a "varying most power- 
ful unfalsified model" (VMPUM) over polynomial Ore algebras such as the Weyl 
algebra or the difference algebra. Mathematically, this amounts to computing 
kernels of module homomorphisms over these algebras. On the one hand, this 
can be achieved using Grobner bases techniques, and on the other, by trans- 
lating the problem to an associated syzygy computation over a commutative 
polynomial ring, thus mimicking ideas of differential algebra. We have also 
studied some structural properties of the resulting models, and we have seen, 
in terms of examples, that models with polynomial coefficients provide a much 
better (and more precise) description of the data than models with constant co- 
efficients. Further future work concerns, for instance, a characterization of the 
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vector space dimension of the VMPUM of several trajectories, thus generalizing 
Theorem 15.21 Let p = [pi, . . . ,p m ] consist of C-linear independent signals. We 
conjecture that dime ( VMPUM (p)) = m. Moreover, it seems possible to us to 
develop VMPUM with polynomial coefficients for data, represented by rational 
and by rational-exponential functions. 
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